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Abstract 

B^B^ flavor oscillations are studied in e'^e~ annihilation data collected with the BABAR detector 
at center-of-mass energies near the T{AS) resonance. One B is reconstructed in a hadronic or 
semileptonic decay mode, and the flavor of the other B in the event is determined with a tagging 
algorithm that exploits the relation between the flavor of the heavy quark and the charges of 
its decay products. Tagging performance is characterized by an efficiency ei and a probability 
for mis-identification, tuj, for each tagging category. We report a determination of the wrong- 
tag probabilities, Wi, and a preliminary result for the time-dependent B^B^ oscillation frequency, 
Amd = 0.512± 0.017± 0.022 ^ps'^ 
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1 Introduction 



The phenomenon of particle-anti-particle mixing in the neutral B meson system was first observed 
almost fifteen years ago [^. The oscillation frequency in B^B^ mixing has been extensively 
studied with both time-integrated and time-dependent techniques Q. In this paper we present 
a preliminary measurement of time-dependent mixing performed at the PEP-II asymmetric e^e~ 
collider, where resonant production of the T{AS) provides a copious source of B^B^ pairs moving 
along the beam axis {z direction) with a Lorentz boost of (3"f = 0.56. The typical separation 
between the two B decay vertices is = P^ctb = 260 /xm, where tb = 1.548 ± 0.032 ps is the B^ 
lifetime The B^B^ mixing probability is a function of Am^, the difference between the mass 
eigenstates Bjj^ and B^, and the time between the B decays, At = Az/P^c: 

Prob{B° W) oc (i _ cos ArudAt). (1) 

In the Standard Model, B^W mixing occurs through second-order weak diagrams involving the 
exchange of up- type quarks, with the top quark contributing the dominant amplitude. A mea- 
surement of Arrid is therefore sensitive to the value of the Cabbibo-Kobayashi-Maskawa matrix Q 
element Vtd- At present the sensitivity to Vtd is not limited by experimental precision on Am^, 
but by other uncertainties in the calculation, in particular the quantity f^Bs, where fs is the 
B^ decay constant, and -B_b is the so-called bag factor, representing the strong interaction matrix 
elements. 

In this analysis, we study the time-dependent probability to observe B^B^, B^B^ and B^B^ 
pairs produced in T(45') decay. We reconstruct one i? in a flavor eigenstate, and use the remaining 
particles from the decay of the other B to identify, or "tag", its flavor. The charges of identified 
leptons and kaons are the primary indicators of the flavor of the tagging B, but other particles 
also carry flavor information that can be identified with a neural network algorithm. The tagging 
algorithm used in this analysis is identical to that employed by BABAR in CP violation studies, in 
which one B is fully reconstructed in a CP eigenstate ||5[. 

Considering the B^B^ system as a whole, one can classify the tagged events as mixed or unmixed 
depending on whether the reconstructed B, referred to as Brec, has the same or the opposite flavor 
as the tagging B, referred to as Btag- If the flavor tagging were perfect, the asymmetry 

/ A .\ -^unmix -^mix /n\ 

«(^^) = AT . +Ar . (2) 

-''^ unmix -"'mix 

plotted as a function of At would describe a cosine function with unit amplitude. However, the 
tagging algorithm incorrectly identifies the tag with a probability w. This mistag rate reduces the 
amplitude of the oscillation by a "dilution factor" 2? = (1 — 2w). When more than one type of 
flavor tag is employed, each will have its own mistag rate, Wi. A simultaneous fit to the mixing 
frequency and its amplitude allows the determination of both Arud and the mistag rates, Wi. 

Neglecting any background contributions, the probability density functions (PDF's) for the 
mixed (— ) and unmixed (+) events can be expressed as the convolution of the oscillatory component 
h±, with a time resolution function 7^(At|d): 

n±{At; r, Amd, V, a) = h±{At; T, Am^, V) ® 7^(A^; a), (3) 



^The symbol refers to the Bd meson; charge conjugate modes are imphed throughout this paper. 
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where a are the parameters of the resolution function and 



h±{At;T, Amd,V) = ^Te-^\^*\ [I ± V cos Ama At] . (4) 

The log-likehhood function is then constructed from the sum of 7i±- over all mixed and unmixed 
events, and over the different tag types, i, each with its own characteristic dilution factor P^: 



tagging 



in£= y: 



J2 lnn+{At-T, Arrid, A, en) + J2 At;r, Am^, A, 

unmixed mixed 



(5) 



The log-likelihood is maximized to extract the dilutions T>i and, simultaneously, the mixing 
parameter Am^. The correlation between these parameters is small, because the rate of mixed 
events at low values of At, where the B^B^ mixing probability is small, is principally governed by 
the mistag rate. Conversely, the sensitivity to Am^ increases at larger values of At; when At is 
approximately twice the B lifetime, half of the neutral Bs will have oscillated. 

Alternatively, the mistag rate can be extracted in a time-independent analysis. Neglecting 
possible background contributions and assuming the B^ec flavor is correctly identified, one can 
express the observed time-integrated fraction of mixed events Xohs ^ function of the B^B^ mixing 
probability Xd ■ 

Xobs = Xd + (1 - 2xd) w, (6) 

where Xd = ^ x^dl + ^d) ^d = Am^^/r. The current world average for Xd is 0.174 it 0.009 Q. 
Because decay time information is available in BABAR, we can improve the statistical precision 
on w by selecting only events that fall into an optimized time interval, |At| < t', where t' is the 
value of I At I above which the integrated number of mixed events equals the integrated number of 
unmixed events. Through the use of an optimized At interval this method achieves nearly the same 
statistical precision on w that is obtained in a full time-dependent likelihood fit. 



2 The BaBAR detector and data set 

The data used in this analysis were collected with the BABAR detector Q at the PEP-II storage 
ring in the period January- June, 2000. The total integrated luminosity of the data set is 8.9 fb~^ 
collected near the T{4S) resonance and 0.8 fb"^ collected 40MeV below the T{4S) resonance (off- 
resonance data). The corresponding number of produced BB pairs is estimated to be about (10.1 ± 
0.4) X 10^. 

The BABAR detector is described in more detail elsewhere Q. For this analysis, the most im- 
portant detector capabilities include charged-particle tracking, vertexing and particle identification. 
Charged particles are detected and their momenta measured by a combination of a central drift 
chamber (DCH) filled with a helium-based gas and a five-layer, double-sided silicon vertex tracker 
(SVT), immersed in a 1.5 T axial field produced by a superconducting magnet. The charged- 
particle momentum resolution is given by (Spt/pt) = (0.0015 pt)^ + (0.005)^, where Pt is m 
GeV/c. The SVT, with typical 10 /xm single-hit resolution, provides vertex information in both the 
transverse plane and in z. The B meson decay vertex resolution is typically 50 /um in z for a fully 
reconstructed B meson and about 100 to 150 //m for the companion (unreconstructed) B meson 
in the event. Beyond the outer radius of the DCH is a detector of internally reflected Cherenkov 
radiation (DIRC) which is used primarily for charged hadron identification. The device consists 
of quartz bars in which Cherenkov light is produced as relativistic charged particles traverse the 
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material. The light is internally reflected, and the Cherenkov rings are measured with an array 
of photomultiplier tubes mounted on the rear of the detector. A CsI(Tl) crystal electromagnetic 
calorimeter (EMC) is used to detect photons, and neutral hadrons, and also for electron identifica- 
tion. The EMC is surrounded by a superconducting coil which provides a 1.5 T magnetic field. The 
Instrumented Flux Return (IFR) consists of multiple layers of resistive plate chambers interleaved 
with the flux return iron and is used in the identification of muons and neutral hadrons. 

3 Particle identification 

Identification of electrons, muons and kaons is an essential ingredient of both B reconstruction 
and flavor tagging. As noted above, the BABAR detector has several systems that contribute to 
particle identiflcation, including the measured dE/dx in the SVT and the drift chamber, Cherenkov 
angle determination in the DIRC, electromagnetic energy measurement in the Csl calorimeter and 
detection of penetrating particles in the IFR. In this section we describe in detail the selections 
that are employed to identify particles by species, both for the purposes of flavor tagging and for 
B reconstruction, where the latter typically employs looser selection criteria. We also include the 
preliminary average values for efficiency and pion misidentification probabilities, as determined 
from data. 

3.1 Electron identification 

For tagging purposes, an electron candidate must be matched to an electromagnetic cluster in the 
Csl calorimeter consisting of at least three crystals, and the ratio of the cluster energy to the track 
momentum E/p must be between 0.88 and 1.3. Cuts @ based on electromagnetic shower shape 
are also used in identifying electron candidates. In addition the electron candidate track is required 
to have a specific ionization (dE/dx) measurement in the drift chamber consistent with that of an 
electron, and, if measured, the Cherenkov angle in the DIRC is required to be consistent with that 
of an ultra-relativistic particle. The electron efficiency and pion misidentification probabilities for 
this selection are about 92% and 0.3%, respectively. 

A somewhat looser electron selection is used in lepton identification for B^ D*~e'^u events. 
In this case, the dE/dx and E/p selections are loosened, while the requirements on electromagnetic 
shower shape and Cherenkov angle are removed. This looser selection has efficiency for electrons 
of approximately 97%, while the pion misidentification probability is about 3%. 

3.2 Muon identification 

Muon identification relies mainly on the number of measured interaction lengths A traversed by the 
track in the IFR iron. At least 2.2A are required, and at higher momenta we require more than Xexp^ 
1, where Xexp is the number of expected interaction lengths as a function of momentum. To reject 
hadronic showers, we make requirements on the number of hit IFR strips that are parametrized as 
a function of the penetration length and the distance of the hits from the extrapolated track. In the 
forward region, which suffers from machine background, extra hit-continuity criteria are applied. 
In addition, if the particle is in an angular region covered by the EMC, the muon candidate must 
have an energy deposit in the calorimeter larger than 50MeV and smaller than 400 MeV. 

These same selection criteria are used for both tagging purposes and muon identification in the 
decay B^ — > D*^ fi^u. The average efficiency and pion misidentification probability are about 75% 
and 2.5%, respectively. 
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3.3 Kaon identification 

Kaon identification employed for flavor tagging requires the ratio of the combined kaon likelihood 
to the combined pion likelihood be larger than 15. The combined likelihoods are obtained by 
multiplying the individual likelihoods from the SVT, DCH and DIRC subsystem information. In 
the SVT and DCH tracking detectors, the likelihoods are based on the dE/dx truncated mean 
measurement compared to the expected value for the K and vr hypotheses, assuming a Gaussian 
distribution. The dB/dx resolution is estimated on a track by track basis given the direction 
and momentum of the track and the number of energy deposition samples. In the DIRC, the 
likelihood is computed by combining the likelihoods computed from the measured Chcrcnkov angle, 
in comparison with the expected Cherenkov angle for a given mass hypothesis, and the Poisson 
probability for the observed number of Cherenkov photons, given the number of expected photons 
for the same hypothesis. DIRC information is not required below 0.7 GeV/c where the DCH dE/dx 
alone provides good KJ-k discrimination. The efficiency of this selection for kaons is about 85% 
and the pion misidentification probability is about 5%. 

In the reconstruction of decays, some channels require a very loose kaon selection to reduce 
backgrounds to acceptable levels. The selection criteria are similar those described above, except 
that the kaon hypothesis is assumed for those subsystems that have no particle identification in- 
formation. This looser selection results in a higher kaon efficiency of about 95% and an acceptable 
pion misidentification probability of 20%. In addition, a loose pion selection is employed in re- 
constructing some B decay modes. This requirement consists of rejecting pion candidates if they 
satisfy the tighter kaon or lepton criteria used for flavor tagging described above. 

4 Time resolution function 

The time difference. At = trec — ttag, between B decays is determined from the measured separation 
Az between the reconstructed B and flavor-tagging decay vertices along the z axis using the known 
boost. At = Az//97c. The At resolution is dominated by the z position resolution for the fiavor- 
tagging B meson vertex. The B^ag production point and three-momentum, with its associated error 
matrix, are derived from the fully reconstructed B^ec candidate three momentum, decay vertex and 
error matrix, and from the knowledge of the average position of the interaction point and the 
T{AS) average boost. These Btag parameters are used as input to a geometrical fit to a single 
vertex, including all other tracks in the event except those used to reconstruct Brec- In order to 
reduce bias due to long-lived particles, all possible F° candidates that can be reconstructed are 
used as input to the fit in place of their daughters. Tracks whose contributions to the are greater 
than 6 are removed from the fit, in an iterative procedure that continues until all remaining tracks 
satisfy this requirement or all tracks are removed. 

The time resolution function can be approximated by a sum of three Gaussian distributions 
with different means and widths 



Fitting the vertex resolution function with simulated events shows that most of the events 
(/i = 1 — /2 — /s ~ 70%) are found in a core Gaussian component, which has a width cri 0.6 ps, 
while the remaining events reside in the tail Gaussian, which has a width (72 ~ 1.8 ps. Tracks from 
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forward-going charm decays included in the reconstruction of the Bfag vertex introduce a smah 
bias, 6i ^ — 0.2ps, for the core Gaussian towards negative values of At. 

A small fraction of events have large values of |At|, due to incorrectly reconstructed vertices. 
This is accounted for in the parameterization of the time resolution function by the third Gaussian 
component, centered at zero with broad fixed width of 8ps, making it almost constant over the 
time interval of the fit. The fraction of events populating this component of the resolution function, 

= /3 is approximately 2 %. 

In the final likelihood fits, we describe the At resolution by introducing two scale factors 5i 
and ^2 that are applied to the event-by-event resolution, a^t, calculated the error on Az provided 
by the vertex fit. We take the width of the core and the tail Gaussian components for each event 
to be (Ti = 5i X a At and a2 = S2 x o"At, respectively. The scale factor Si and the bias 61 of the 
core Gaussian are free parameters in the fit. The scale factor ^2 and the fraction of events in the 
core Gaussian /i are fixed to the values estimated from Monte Carlo simulation by a fit to the pull 
distribution (^2 = 2.1 and /i = 0.75). The bias of the tail Gaussian, 62, is fixed at Ops. The three 
free parameters in the likelihood fit are 

a = {Si, 61, U} . (8) 

We observe no significant differences in simulated events between resolution function parameters 
obtained from samples involving different decay modes for Brec- This is expected, because At 
resolution is dominated by the precision on the Bfag, rather than the Brec vertex. Likewise, the 
differences in the resolution function parameters for the different tagging categories are also small. 
Therefore, we use a single set of resolution function parameters a for all decay modes. Table |l| 
shows the values for the vertex parameters obtained in data from a fit to the hadronic B^ sample, 
described below. 

Table 1: Parameters of the resolution function determined from the sample of events with hadronic 
fully reconstructed B candidates. 



parameter 


value 


61 (ps) 


-0.20 ±0.07 


from fit 


Si 


1.33 ±0.13 


from fit 


U (%) 


1.6 ±0.6 


from fit 


/i (%) 


75 


fixed 


52 (ps) 





fixed 


S2 


2.1 


fixed 



5 Flavor tagging 

After the daughter tracks of the reconstructed B are removed, the remaining tracks are analyzed 
to determine the flavor of the Btag, and this ensemble is assigned a tag flavor, either B^ or B^. 

To illustrate the tagging discriminating power of each tagging category, we use as a figure 
of merit the effective tagging efficiency = x (1 — 2u)j)^, where is the fraction of events 
associated to the tagging category i and Wi is the mistag fraction, the probability of incorrectly 
assigning the opposite tag to an event of this category. 
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We use four different types of flavor tag, or tagging categories, in this analysis. Two of these 
tagging categories rely upon the presence of a prompt lepton or one or several charged kaons in the 
event. The remaining two categories, called neural network categories, are based upon the output 
values of a neural network algorithm applied to all the events that have not already been assigned 
to one of the Lepton or Kaon tagging categories. 

5.1 Lepton and kaon tags 

The Lepton and Kaon tagging categories use the correlation between the charge of a primary lepton 
from a semileptonic decay or the charge of a kaon, and the flavor of the decaying b quark. For 
the Lepton category we use both electrons and muons. A minimum ccntcr-of-mass momentum 
requirement on the lepton is applied to reduce the contamination from softer opposite-sign leptons 
coming from charm semileptonic decays. There are no such discriminating kinematic quantities 
to reduce the contamination of opposite-sign kaons, so the optimization relies principally on the 
balance between the kaon identification efficiency and the purity of the kaon sample. 

A lepton tag is defined by taking the charge of the fastest identified electron or muon with a 
center of mass momentum greater than 1.1 GeV/c. A kaon tag is defined by taking the sum of the 
charges of all identified kaons. If both an electron and a muon are identified, the electron tag takes 
precedence. If the event has a lepton tag and there is no conflicting kaon tag, the event is assigned 
to the Lepton category. If the event has no lepton tags but has a non-zero kaon tag, the event is 
assigned to the Kaon category. If the event has both lepton and kaon tags but they conflict, the 
event is not assigned to either the Lepton or the Kaon category. 

5.2 Neural network tags 

The use of a second tagging algorithm is motivated by the potential flavor-tagging power carried 
by non-identified leptons and kaons, softer leptons from charm semileptonic decays, soft pions from 
D* decays, and more generally, by the momentum spectrum of charged particles from B meson 
decays. The best way to exploit the information contained in a set of correlated quantities is to 
use multivariate methods, such as neural networks, rather than to apply selection cuts. 

We design five different neural networks, called subnets, each with a specific goal. Four of the 
subnets are track-based: the L and LS subnets are sensitive to the presence of primary and cascade 
leptons respectively, the K subnet to that of charged kaons and the Sof tPi subnet to that of soft 
pions from D* decays. In addition the Q subnet exploits the charge of high-momentum particles in 
the event. 

The L and LS subnets share a set of discriminating kinematic variables in addition to the boolean 
outputs of the standard lepton identification algorithms. The variables are designed to more fully 
characterize the kinematic features of a true primary lepton from a semileptonic B decay. Taking 
all charged tracks in the event, excluding that under consideration as the lepton, together with all 
neutral clusters we form the recoiling system X. Assuming that Btag is produced at rest in the 
center-of-mass frame, and that the excluded track is a primary lepton from the Btag semileptonic 
decay, one can calculate the four-momentum of the neutrino pi, and the four-momentum pw of 
the virtual W boson in the semileptonic decay. The variables available are the center-of-mass 
momentum of the lepton candidate track, the mass and energy of the recoil system X, the cosine 
of the angle between the excluded track momentum and the direction of the neutrino, the cosine of 
the angle between the W momentum direction and the direction of the closest charged or neutral 
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candidate with energy larger than 50 MeV, and the total energy flow in the hemisphere around the 
W momentum direction. 

The K subnet uses the particle momentum in the laboratory frame, together with the three 
relative hkelihoods Lx/iL-n + Lk) for the SVT, the DCH and the DIRC. 

The SoftPi subnet uses the track with the lowest momentum from the Btag decay and the 
angle it makes with respect to the thrust axis, calculated using all other charged tracks. 

The variables of the Q subnet are the momentum of the highest momentum track in the Btag 
system and the sum of momentum- weighted charges normalized by the momentum sum of all tracks 
with impact parameter less than 1 mm and all neutral clusters with energy greater than 50 MeV. 

The variables used as input to the neural network tagger are the highest values of the L, LS 
and SoftPi subnet outputs each multiplied by the charge of the corresponding tagging track, the 
highest and the second-highest values of the K subnet output again multiplied by the charge of the 
corresponding tagging tracks, and the output of the Q subnet. 

The output from the full neural network tagger, xmt^ can be mapped onto the interval [—1,1]. 
The assigned flavor tag is if xnt is negative, and B^ otherwise. Events with Ixat^I > 0.5 are 
assigned to the NTl tagging category and events with 0.2 < \xnt\ < 0.5 to the NT2 tagging category. 
Events with |xArr| < 0.2 have very little tagging power and are rejected. 

6 Event selection and sample composition 

B^ mesons are reconstructed in the hadronic decay modes B^ D^*^~'it~^, D^*'^~ , D^*^~af, 
J/tp K*^ and the semileptonic decay mode B^ — > D*~£~^i'. All final state particles, with the excep- 
tion of the neutrino in the semileptonic decay, are reconstructed. 

Charged tracks measured by the drift chamber and/or SVT are required to originate within 
1.5 cm in xy and 10 cm in z of the nominal beamspot. Electromagnetic clusters in the calorimeter, 
that are unassociated with charged tracks and used to reconstruct vr^ candidates, are required to 
have an energy greater than 30 MeV and a shower shape consistent with a photon interaction. 

Neutral pion candidates are formed from pairs of electromagnetic clusters assumed to be pho- 
tons. The invariant mass of the photon pair is required to be within ±20 MeV/ (2.5 a ) of the 
nominal mass, with a minimum summed energy of 200 MeV. Selected candidates are kinemati- 
cally fitted with a mass constraint. 

K'^ candidates are reconstructed in the tt'^tt^ mode, with an invariant mass, computed at the 
vertex of the two tracks, between 462 and 534 MeV/ c^. The of the topological vertex fit is 
required to have a probability greater than 0.1%. The opening angle between the flight direction 
and the momentum vector for the Kg candidate is required to be smaller than 200 mr. Finally, the 
transverse flight distance from the primary vertex in the event, rxy, is required to be greater than 
2 mm. 

In order to reject "jet-like" events from e"^e^ qq (continuum) background, we require that the 
normalized second Fox- Wolfram moment (i?2 = H2/Hq), calculated with both charged tracks 
and neutral clusters, be less than 0.5 (0.45) in hadronic (semileptonic) decay modes. 

6.1 Hadronic decays 

B^ mesons are reconstructed in the hadronic modes B'^ D^*^~Tr~^ , B^ D^*^~ , B^ D^*^~af 
and B^ —>■ J/ipK*^. A variety of , D~ , and D*~ modes are used to achieve reasonable efficiency 
despite the typically small branching fractions for any given B decay channel. 
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6.1.1 Event selection 

We select D^, D~ and D*~ mesons with the foUowing criteria. candidates are identified in the 
decays channels K^tt'^, K^7:~7:^, it~^ it~ 'k~ and K^'k^-k~ . D~ candidates are selected using 
the K^t:^t:^ and K^tt^ modes. The D*^ candidates are found using the decay D*~ — > D^7r~. 
Charged and neutral kaons are required to have a momentum greater than 200 MeV/ c. The same 
criterion was applied to the pion in — > D^*^~Tr~^, — > D^*^~ decay. For the decay modes 
B^ — > D^*^~af, the pions are required to have momentum larger that 150MeV/c . and D~ 
candidates are required to lie within itScr, calculated on an event-per-event basis, of the nominal 
masses. Pull distributions for the and meson masses have been measured in data and were 
found to have rms in the range of 1.1-1.2 when fitted to a Gaussian form. For — > K^it~tt^, we 
only reconstruct the dominant resonant mode — > K~^p~, followed by p~ — > 7r~7r'^. The tt~tt^ 
mass is required to lie within ±150 MeV/ of the nominal p mass and the angle between the vr" 
and in the p rest frame, O'^o^, must satisfy | cos^IjOtt I> Finally, all and D~ candidates 
are required to have a momentum greater than 1.3GeV/c in the T(4S) frame and a probability 
of the topological vertex fit greater than 0.1%. A mass constrained fit is applied to candidates 
satisfying aforementioned requirements and is used in the subsequent reconstruction chain. 

We form D*~^ candidates by combining a with a pion which has momentum greater than 
70MeV/c. The soft pion is constrained to originate from the beamspot when the D*~ vertex is 
computed |^]. To account for the small energy release in the decay T(4S') — > BB (resulting in a 
small transverse flight of the B candidates), the assumed vertical size of the beam spot is inflated 
to 40 /im. Monte Carlo simulation was used to verify that this inflation does not introduce any 
signiflcant bias in the selection or in the At measurement. D*^ candidates are then required to 
have Am = m{D^-K~) — m{D^) within ±1.1 MeV/c^ of the nominal value for — > K~^'k~it^ mode 
and ±0.8 MeV/ for all the other modes. This corresponds to about ±2.5(T, where the resolution is 
estimated by taking a weighted average of the core and broad Gaussian components of the observed 
Am distributions. 

B^ candidates are formed by combining a D*^ , D~ or J/-0 candidate with a tt^, p^, or K*^. 
For B^ D*~p^, the tt^ from the p'^ decay is required to have an energy greater than 300 MeV. 
For B^ D*^ai, the af is reconstructed by combining three charged pions, with invariant mass 
in the range 1.0 to 1.6 GeV/c^. In addition, the probability of the vertex fit of the af candidate 
is required to be greater than 0.1%. 

As described in Section |3|, kaon identification is used to reject background. For most B^ modes, 
no identification or only a loose selection is enough to achieve signal purities of 90%. In the B^ — > 
D~ai mode reconstruction, a tight kaon identification is required to alleviate large combinatorial 
backgrounds. 

In order to suppress continuum background, in addition to the R2 requirement, we calculate 
the thrust angle, Ofh, defined as the angle between the thrust axis of the particles which form the 
reconstructed B^ candidate and the thrust axis of the remaining tracks and unmatched clusters 
in the event, computed in the T(4S) frame. The two thrust axes are almost uncorrelated in BB 
events, because the B^ mesons are almost at rest in the T(4S') rest frame. In continuum events, 
which are more jet-like, the two thrust axes tend to have small opening angles. For final states with 
a D* and 2 (3) pions we require | cos6th\ < 0.9 (0.8) for the — > K~^7r~ and K~^7:~7r^ modes and 
0.8 (0.7) for the — > K~^7r~^7:~'K~ and i^g7r"'"7r~. No such requirement is applied in the case of 
the B^ D*~7r~^ mode. In modes which contain a D and a single (2,3) pion(s) in the final state, 
we require |cos6'(ft| < 0.9 (0.8, 0.7). 

The B° J/'4)K*^ selection is the same as described in Ref. The J/if: e'^e candi- 
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dates are formed from pairs of oppositely-charged tracks assumed to be electrons, as explained in 
Section ^ At least one of the decay products must be positively identified as an electron or, if 
outside the acceptance of calorimeter, must be consistent with the electron hypothesis from the 
dE/dx measurement in the drift chamber. If both tracks are in the calorimeter acceptance and 
have a value of E/p larger than 0.5, an algorithm for the recovery of Bremsstrahlung photons Q 
is used. The J/ijj — > candidates are formed from oppositely-charged tracks identified as 

muons, as described in Section |^. At least one of the decay products must be positively identified 
as a muon, and the other, if in the acceptance of the calorimeter, must be consistent with a min- 
imum ionizing particle. We retain only those J/ip candidates with an invariant mass in the range 
2.95 < m{J/i;) < S.UGeV/c^ for the e+e" mode and 3.06 < m(J/V') < S.UGeV/c^ for the n+fi' 
mode. We reject events with cos 9th greater than 0.9. 

candidates are characterized by a pair of nearly uncorrelated kinematic variables, the dif- 
ference between the energy of the B^ candidate and the beam energy in the T{4S) center-of-mass 
frame, AE, and the beam energy substituted mass, m-Es 0- In the (mES)^-E') plane the sig- 
nal region is defined as 5.270 < rriEs < 5.290 GeV/c^ and \A.E\ < 3a ae- The sideband region 
is defined as 5.2 < mEs < 5.26GeV/c^ and \AE\ < Sa^E- The value of a^E, the AE resolu- 
tion, is mode-dependent and varies between 19 to 40 MeV. When multiple B^ candidates (with 
mES > 5.2 MeV/ c^) are found in the same event, the candidate with the smallest value of \AE\ is 
selected. 

6.2 Sample composition 

We use the rriES sideband region to measure the combinatorial background fraction, to charac- 
terize the background At distribution, and to measure the fraction of mixed events contributed 
from background under the B^ signal peak. It is therefore useful to check that the background 
composition in the mES sideband region is similar to that in the signal region. The validity of this 
background estimation procedure was checked on 2.0 fb~^ of simulated data. We also compared 
the mES sideband shape in data with the Monte Carlo and found good agreement. 

The B^ signal yield and sample purity extracted from fits to the uies distribution are summa- 
rized in Table |2|. The net B^ signal sample, before applying any decay vertex requirements, consists 
of 2577± 59 signal candidates with a purity of about 86%. 

Table 2: Two-body hadronic B^ decay candidate yields and signal purities from the fit to the rriES 
distribution. Signal purities are estimated for uies > 5.27GeV/c^. 



Decay mode 


Number of B^ candidates 


S/(S+B) [%] 


B'> D*-TT+ 


622 ± 27 


90 


B° D*-p+ 


419 ± 25 


84 


B° D*-at 


239 ± 19 


79 


BO ^ D-TT+ 


630 ± 26 


90 


B^ D-p+ 


315 ± 20 


84 


5° ^ D-al 


225 ± 20 


74 


B° J/iPK*'^ 


194 ± 15 


90 


Total 


2577 ± 59 


86 



In order to select B" candidates with well understood B decay vertex error, we require the 
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convergence of the topological fit for decay vertex. candidates with |Az| > 3.0 mm or with 
CAz > 400 fim are removed. 

The efficiency for each tagging category is calculated as the number of signal events for each 
tag, divided by the total number of signal events after vertex cuts are imposed. The tagging 
efficiency and signal purity for the individual tagging categories in data are extracted from fits to 
the rriEs distributions shown in Fig. |l| and are listed in Table ^. The distributions are fitted with 
a Gaussian distribution for the signal and the ARGUS background function for the background 
parametrization All fits have good confidence levels. 

Table 3: Tagging efficiencies for hadronic B^ decays and signal purities in data separately for the 
four tagging categories. Signal purities are estimated for uies > 5.27GeV/c^. 



Tagging Category 


Efficiency [%] 


B candidates 


S/(S+B) [%] 


Lepton 


10.5 ± 0.6 


260 ± 17 


95 


Kaon 


36.7 ± 1.0 


918 ± 34 


86 


NTl 


12.0 ± 0.7 


305 ± 19 


89 


NT2 


16.4 ± 0.7 


405 ± 24 


81 


All tags 


75.6 ± 0.9 


1886 ± 49 


87 



6.3 Semileptonic B decays 

The semileptonic decay B^ — > D*~i~^i', with a measured branching fraction of 4.6 it 0.27% is 
a copious source of B^ mesons. We reconstruct the D*~ through its decay to D^tt", and use the 
three decay modes K~^7r~, K^t:~^'k~'k~ and K^tt'-k^. 

6.3.1 Event selection 

We reconstruct candidates in the three modes listed above. All reconstructed candidates are 
required to have an invariant mass within ±2.5(7 of the nominal mass. The topological vertex 
fit is required to have a probability greater than 1%; for the — > K^it~t:^ mode this vertex 
fit includes a kinematic constraint on the vr'' mass. There are no additional requirements for 
K~^TT~ . For K^TT^TT^TT^ and — > K^ir^ir^ we require loose K and vr particle identification 

as described in Section |3|, and a minimum vr*^ momentum of 200 MeV/ c. In addition, the K and tt 
candidates are required to have momenta greater than 200 and 150 MeV/ c, respectively, for the mode 

— > K~^7:~^TT~7r~ . The decay K~^7r~Tr^ occurs mostly through resonant substructures. The 

p and K* resonances dominate and we use the measured Dalitz weights to construct an event-by- 
event probability and select events using this quantity to suppress combinatorial background]^. 

candidates satisfying the above requirements are combined with all charged tracks, having a 
minimum transverse momentum of 50MeV/c and charge opposite to that of the candidate kaon, to 
form D* candidates. The mass difference m{D*~) — m{D^) is calculated and required to lie within 
2.5(7 of the nominal value. 

Finally, D* candidates are combined with electron or muon candidates, with momentum greater 
than 1.2GeV/c and satisfying the lepton identification requirements described in Section ^. The 

^ The Dalitz weights are used in the time-integrated analysis. The time-dependent analysis uses the — > K'^ p~ 
sample only. 
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Figure 1: tties distribution for each tagging category (Lepton, Kaon, NTl and NT2) for all the 
hadronic modes. 
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D* - lepton topological vertex fit is required to have a probability greater than 1%, and the 
lepton and D* candidates must have opposite charge. The D* and lepton tend to be back-to-back 
in the rest frame, so we require cos 9{D* — I) < where 9{D* — I) is the angle between the D* 
and the lepton in the center-of-mass frame. 

The neutrino cannot be reconstructed, but we require that the candidate four- momenta be 
consistent with a missing particle of zero mass: 

{PB - PD* -Pif =pI = 0. (9) 

Solving this equation in the T(4S') frame, we obtain a constraint on the angle between the and 
the D* — I system: 

„,,,B.„.,).M±MiZ^f£^. (10) 

2|pbIIpd.iI ^ ' 

Only the mass and scalar momentum of the initial state S^, both of which are known, is required 
to calculate the angle 9{B, D*l) and not the flight direction. In B^ D*lv decay, the cosine of this 
angle must lie in the region (— 1,+1). Allowing for detector resolution effects in the reconstructed 
momenta and angles, we require | cos6{B — D*l)\ < 1.1. 

After applying these criteria, we obtain a sample of 7517ibl04 B D*lv events: 3101 ±64 in the 
Z)° K+TT- mode, 1986 ±51 in the Z)° K+-k-'k^ mode, and 2430 ±56 in the Z)° K+tt+tt-tt- 
mode. 



6.3.2 Sample composition 

The events are flavor tagged as described in Section |5|. The D* — mass difference distributions 
for each tagging category are shown in Fig. |2[ The backgrounds are evaluated separately for each 
tag category. Backgrounds are larger for the semileptonic modes than for the hadronic modes and 
originate from a variety of sources. Each source of background is evaluated using a control sample 
that is taken from data whenever possible. The background control samples are used to measure the 
background fractions, to characterize the background At distribution, and to measure the fraction 
of mixed events contributed by each source of background. 

We divide the backgrounds to B^ D*lv into three types: events with an incorrectly recon- 
structed D* ("combinatorial" background), events in which a true D* is combined with an incorrect 
lepton candidate ( "wrong- lepton" background), and events in which a correctly identified -D*-lepton 
pair originates from semileptonic B^ decay, for example from the decay B^ D**^l~^u, followed 
by D**^ D*~mT^ ( "i3+ background"). The symbol D**^ refers to an admixture of orbitally or 
radially excited charmed meson resonances which decay via strong interaction into a D*~ . Events 
of the type B^ — > D**~l'^v are considered as signal. 



Combinatorial background 

The fraction of combinatorial background, due to falsely reconstructed D* candidates, is esti- 
mated by fitting the Am(Z)* — D^) distributions. We use a Gaussian distribution to characterize 
the signal and a threshold function with a sharp rise followed by an exponential tailoff to char- 
acterize the background shape, the signal region is defined to lie within ±2.5(T of the peak in 
Am{D* - D°), while events in the the sideband region 0.150 < Am{D* - I?°) < 0.160 GeV/c^ 
provide the combinatorial background control sample. 
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Figure 2: D* — mass difference distribution for each tagging category (Lepton, Kaon, NTl and 
NT2) for the 0*11^ sample. 
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Wrong-lepton background 

There are four potential sources of background events in which a real D* is combined with a 
wrong lepton. We consider each in turn. 

First, there are events in which the lepton candidate is not a real lepton, but is misidentified as 
such. This "fake lepton" background is estimated from data by selecting events in which a track 
candidate that has failed very loose lepton criteria is substituted for the lepton candidate. These 
events are weighted with the lepton mis-identification probabilities measured in data to estimate the 
fraction of fake lepton background, after subtraction of the combinatorial background contribution. 

The second type of wrong lepton events is due to a real D* from one B combined with a 
real lepton from the other B. These "uncorrelated lepton" backgrounds are due to B^B^ events 
where mixing has occurred, yielding the right D*l sign combination, or events in which a lepton 
originating from a secondary charm decay is combined with a D* from the other B. The D* 
and lepton directions are uncorrelated in these events, prior to application of the cos 9{D*l) and 
cos 6{B — D*l) selection criteria. 

To estimate this background from data, the lepton momentum vector in the T{AS) center-of- 
mass frame is parity-inverted, prior to the calculation of the 0{D* — I) and 6{B — D*l). Events with 
a correlated back-to-back D* — I pair fail these criteria after the lepton momentum is flipped, but 
events with a randomly correlated D* — I pair pass with approximately the same efficiency as in 
the original sample. After removing the remaining signal contribution, and correcting for residual 
combinatorial background, we estimate the uncorrelated lepton background from the fiipped-lepton 
control sample. 

Third, there are decays of the type B^ D*DX, where the D decays semileptonically produc- 
ing a non-primary lepton. The momentum requirement on the lepton rejects most of these "cascade 
lepton" events; the remaining fraction, estimated from Monte Carlo simulation, is less than 1% and 
has been neglected. 

Finally, cc events can produce a real D* and a real lepton in a back-to-back configuration. The 
cc background fraction has been estimated using combinatorial-subtracted off-resonance data. We 
use the off-resonance data as a control sample to characterize the At distribution and the mixed 
event contribution from this background source. 

S+ background 

In addition to the well-studied semileptonic decays B — > Dlv and B D*lu, a significant 
fraction of B semileptonic decays involve additional final state particles, either produced through 
D** resonances or non-resonantly. These processes of the type B — > D*{mT)lv contribute to both 
the neutral and charged semileptonic B decays. 

For the purposes of this analysis, neutral semileptonic decays B^ D*[n'K)lv that pass our 
event selection criteria are considered to be part of the signal. They contribute equally to the 
measurement of Am.fi, and the additional low- momentum pion does not affect the tagging algorithm. 

However, the charged B decays of the type B~ — > D*^ {mT)l~ v are considered as background 
for this analysis. They do not oscillate and must be corrected for in extracting Am^; their mistag 
rate may differ from that of B^ decays as well. 

We assume a lifetime and At resolution for charged B events that is the same as that for the 
B^ signal events. The mistag rate is estimated directly from data using fully reconstructed B^ 
decays. It is difficult to accurately estimate the fraction of this background, since it cannot be 
cleanly separating from the signal. To estimate the fraction of events due to charged B decay. 
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Table 4: Sample composition and yields in tagged B" D*lv events 





Lepton 


Kaon 


NTl 


NT2 


Combinatorial 
Fake Lepton 
Uncorr. Lepton 

cc 

Signal Fraction 


0.063±0.008 
0.032±0.007 
0.015±0.016 
O.OOOibO.007 
0.062±0.040 
0.827±0.040 


0.161±0.007 
0.035±0.007 
0.035±0.012 
0.026±0.010 
0.052±0.031 
0.691±0.031 


0.128±0.011 
0.031±0.007 
0.024±0.013 
0.025±0.017 
0.055±0.038 
0.737±0.038 


0.158±0.009 
0.030±0.007 
0.028±0.013 
0.056±0.021 
0.051±0.034 
0.677±0.034 


Tagged Event Yield 
Tag Efficiency 


863 ±32 
0.121 ±0.007 


2804 ±63 
0.368 ±0.017 


850 ±34 
0.114 ±0.006 


1318 ±43 
0.169 ±0.009 



we rely on an estimate of their production rate, combined with a Monte Carlo study of the event 
selection efficiency for these modes. 

The inclusive branching fraction, B{B^ D*^ {mT)l'^v) = 1.25±0.16% is taken from LEP 



measurements |1C]. We use Monte Carlo simulation to study the efficiency of the event selection 
criteria for events of this type, averaging over several resonant D** and non-resonant states. The 
average efficiency is found to be 4.25± 3.0%, where a conservative systematic error has been assigned 
due to the lack of knowledge of the relative decay fractions for the possible modes. The product of 
the branching ratio times the event selection efficiency, B x e{B~^)^ and the corresponding product 
B X e{B^) for the neutral B signal events are computed, where decays of the type B^ — > D*^{mT)l~i' 
are included as part of the signal. By this method, we find: 



Background summary 

The background contributions in i? ^ D*lv events, averaged over decay modes, are sum- 
marized by tag category in Table ^. 

7 Likelihood fit method 

In the presence of backgrounds, the probability distribution functions l-i± of Eq. ^ must be extended 
to include a term for each significant background source. The background parameterizations are 
allowed to differ for each tag category. Each event is identified as being either mixed (— ) or unmixed 
(±) and as belonging to a particular tag category, i. Thus a distribution must be specified for each 
possibility (±/— 

n±,^ = fsig,in±{M-T,^mD,Vi,ai)+ ff^,i>3±,id^t;b±^i^p) (12) 

/3=bkgd 

where the background PDFs, B±^i^p, provide an empirical description the At distribution of the 
background events in the sample. The fraction of background events for each source and tagging 
category is given by fp^i, while ?)±,i,/3 are parameters used to characterize each source of background 
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by tagging category for mixed and unmixed events. The signal fraction in each tag category is given 

by 

fsig,i = 1 - XI fP'i (^^) 

The distributions are normahzed so that for each i and P 

/oo 
dAtiB+,,,p + B.,,,p) = 1. (14) 

7.1 Background parameterization 

Backgrounds stem from many different sources; we use control samples, derived whenever possible 
from the data itself, to characterize the background time dependence and dilution. These control 
samples were described in Section ^. We use an empirical description for the time dependence of 
the backgrounds in the likelihood fit allowing for three time components for each background, each 
with its own dilution factor D and a common resolution function TZ : 

• Zero lifetime component: -B±,i = (l=bD^)(5(t) (8) TZ 

• Non-zero lifetime component, no mixing: -B±,2 = ^(li^2) e^"'"^'^*'-' (^TZ 

• Non-zero lifetime component, with mixing: B± 3 = ^^e''~'"^l^*l^(libD3 cos(Am3At)) 7?, 

A likelihood fit to the background control samples is used to determine how much of each time 
component is present and to fit for the apparent lifetimes, resolution, mixing and dilutions that best 
describe the background sample. This approach allows for more fit parameters than are absolutely 
necessary. The goal is to determine the background shapes as well as possible in an empirical sense. 



7.2 Likelihood fit results 
7.2.1 Hadronic decay modes 

We extract Am^ and the mistag rates by fitting the At distributions of the selected B candidate 
events with ttt-es > 5.2 GeV/c^ with the likelihood function described above. The probability that a 
B candidate is a signal or a background event is determined from a fit to the mES distribution. We 
describe the rriEs shape with a single Gaussian S'(mEs) for the signal and the ARGUS function Q 
-S(mEs) for the background. Based on this fit, the event-by-event signal probability is determined 
from: 

= Sirn'sHBirn^s) ^''^ 
The contribution of each event to the fitted signal parameters corresponds to this signal probability. 

The At distributions of the combinatorial background are described with a zero lifetime com- 
ponent and a non-oscillatory component with non-zero lifetime. We fit for separate resolution 
function parameters for the signal and the background in order to minimize correlations between 
the background parameters and the signal parameters. 

We use the data sample described in Section ^ consisting of « 1900 fully-reconstructed and 
tagged B^ candidates. The At distributions of those candidates, overlaid with the likelihood fit 
results, are shown in Fig. ^ for the candidates with tit-es > 5.27GeV/c^. The At distributions 
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for the background candidates in data from the m-Es sideband (rn-Es < 5.27GeV/c^) are shown 
m Fig. |. In Fig. | the mixing asymmetry of Eq. ^ is plotted; the primordial cosine function is 
clearly visible. The results from the likelihood fit to data are summarized in Table ^. The tagging 
separation Q = etag{^ — 2w)'^ is calculated from the mistag rate and the efficiency quoted in Table |3|. 
Summing over all tagging categories, we measure a combined effective tagging efficiency Q ~ 28%. 

Table 5: Results from the likelihood fit to the At distributions of the hadronic and semileptonic 
B decays. Amrf and the mistag rates include small corrections corresponding to the difference 
between the generated and reconstructed values in simulated signal events. The summed Q over 
all tagging catagories is 0.285 (0.283) for hadronic (semileptonic) decay modes. 



Parameter 


Hadronic 




Semileptonic 




Fit Value Q = 


e(l-2u')2 


Fit Value Q 


= e{l-2wf 


Am^ [n ps"^] 


0.516 ± 0.031 




0.508 ± 0.020 




ii;(Lepton) 


0.116 ± 0.032 


0.062 


0.084 ± 0.020 


0.071 


ti;(Kaon) 


0.196 ± 0.021 


0.136 


0.199 ± 0.016 


0.133 


tt;(NTl) 


0.135 ± 0.035 


0.064 


0.210 ± 0.028 


0.066 


■u;(NT2) 


0.314 ± 0.037 


0.023 


0.361 ± 0.025 


0.013 


SCalCcore, sig 


1.33 ± 0.13 




1.32 ± 0.07 




'^core, sig [ps] 


-0.20 ± 0.07 




-0.25 ± 0.04 




/outlier 


0.016 ± 0.006 




0.000 ± 0.002 





As a validation of these results we have carried out identical analysis procedures on simulated 
Monte Carlo events generated with a detailed detector simulation, processed through the event 
reconstruction chain in the same manner as the data. In the signal simulation, the fitted values 
of the B^B^ oscillation frequency, Am^^ = 0.451ib 0.011 ?ips^^ (hadronic decays) and Am^ = 
0.456ib 0.009 ?ips~^ (semileptonic decays), are consistent with the value of 0.464 ?ips~^ used for 
Monte Carlo generation. The fitted tagging dilutions and mistag rates are in good agreement with 
the values obtained from Monte Carlo truth information, confirming an unbiased measurement of 
those parameters. We apply the observed differences as a correction to the measured values in 
data, although all parameters as determined from the simulated sample are consistent with the 
generated values. 

7.2.2 Semileptonic decays 

The event selection for the B D*lv sample is summarized in Section ^ The yield of tagged 
B — > D*lv events and the signal purity are described in Section |6[ 

The measurement of B^B^ mixing and the extraction of Am^ and the mistag fractions with 
B — > decays proceeds in two steps. First, we fit the background control samples described in 
Section |6| and determine their parameters. 

Second, we fit the signal events, fixing the background fractions to the values summarized in 
Table § and the At resolutions and dilution parameters to the values obtained from fits to the 
control samples. The results of the fit and the calculated tagging performance, Q = etag (1 - 2wf, 
are summarized in Table ^. The At distribution of the B D*lv candidates, overlaid with the 
likelihood fit result, is shown in Fig. |6| 
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T 1 1 1 1 1 1 r 




A? (ps) 

Figure 3: At distributions in data for the signal hadronic B sample with tties > 5.27GeV/c^. 
The fitted At shapes for the selected candidates and for the fraction of background candidates are 
overlaid. The confidence level for this projection of the fit result is calculated from the binned At 
distributions using a Poisson-x^ technique to be 13% . 
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150 




-8 -4 4 8 

A? (ps) 

Figure 4: At distributions in data for the background candidates with mss < 5.27 GeV/ separately 
for unmixed and mixed candidates. The fitted At shapes for those candidates are overlaid. The 
confidence level for this projection of the fit result is calculated from the binned At distributions 
using a Poisson-x^ technique to be 21% . 
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Figure 5: Time-dependent asymmetry a{At) between unmixed and mixed events for hadronic B 
candidates with m-ES > 5.27GeV/c^. 
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600 




A? (ps) 

Figure 6: At distributions in data for the selected B D*lu candidates for unmixed and mixed 
candidates. The fitted At shapes for those candidates and for the fraction of background candidates 
are overlaid. The confidence level for this projection of the fit result is calculated from the binned 
At distributions using a Poisson-x^ technique to be 28%. 
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2 4 6 8 

At(ps) 

Figure 7: Time-dependent asymmetry between unmixed and mixed events for B D*lu candi- 
dates. 
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7.3 Systematic error estimation 

The systematic errors for the hadronic and semileptonic B samples are summarized in Tables |^ 
and 1^ and can be grouped into three categories: 

At reconstruction 

We determine the level of any potential systematic bias in Am^ and the mistag rates due to 
the At resolution function by independently varying the scale factor and the mean of the wide 
Gaussian and the fraction of events in the wide Gaussian. The variations correspond to a change 
in the RMS of the total resolution function by one standard deviation as measured in data. The 
systematic uncertainty due to At outliers is estimated from the variation of the fitted parameters 
with the fraction of outliers and the RMS of their At distribution in a sample of toy Monte Carlo 
events. An error in the boost of the T(4S) system or in the z scale of the detector can bias the 
Anid measurement because those parameters are used to reconstruct a decay length difference Az 
and to convert it to the decay time difference At. In the likelihood fit, we fix the lifetime to the 
PDG value ||3| and its uncertainty leads to a systematic error. 

Background parameters 

The signal probability assigned to each candidate in the hadronic B sample has a statistical 
uncertainty, and these statistical uncertainties lead to systematic uncertainties in Am^ and the 
mistag rates. We estimate these uncertainties by varying the width and height of the fitted peak 
in rriEs by one standard deviation. 

In the semileptonic sample we vary the average background fractions by the statistical uncer- 
tainties derived from the control samples. The uncertainty in the fraction of D*ll.u candidates with 
a fake lepton includes the uncertainty in the lepton identification rates. To estimate the sensitivity 
of Anid and the mistag rates to the At description of the combinatorial background, we repeat the 
likelihood fit with an additional oscillatory term in the PDF. In the fit to the signal sample we 
vary the background dilutions obtained from the control samples by one standard deviation. We 
study the sensitivity to the resolution function of the combinatorial backgrounds by allowing an 
additional scale factor to account for possible tails in the At distribution. 

Check with simulated events 

Candidate selection criteria can cause systematic biases in the measurement of Am^ and the 
mistag rates. These biases are estimated with fully simulated events and are found to be consistent 
with zero within their statistical uncertainty. Nevertheless, we correct for the actual differences 
between the generated and reconstructed values and include the statistical uncertainties in the 
measured parameters from the simulated events as a contribution to the systematic uncertainties 
in the result. 

8 Time-integrated method 

8.1 Description of the method 

As previously described in Section ||, a time-integrated method to measure the mistag fractions 
in data provides a simple and robust check of the likelihood fit method presented in Section |^. 
The statistical precision of this method is enhanced by restricting the sample to events in a single 
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Table 6: Systematic uncertainties for Am^ and the mistag rates measured with hadronic B decays 
for the hkehhood fit. 



Source 


Arud 


Lepton 


Kaon 


NTl 


NT2 




[h ps-i] 










At Resolution 


0.011 


0.004 


0.004 


0.004 


0.004 


Background At 


0.002 


0.002 


0.002 


0.002 


0.002 


Background Resolution 


0.002 


0.002 


0.002 


0.002 


0.002 


Background Fractions 


0.004 


0.004 


0.002 


0.006 


0.004 


lifetime 


0.005 


0.001 


0.001 


0.001 


0.001 


z scale 


0.005 










z boost 


0.003 










Monte Carlo Correction 


+0.013 


-0.001 


0.000 


-0.010 


-0.015 




± 0.011 


± 0.011 


± 0.008 


± 0.015 


± 0.014 


Total Systematic Error 


0.018 


0.013 


0.010 


0.017 


0.015 


Statistical Error 


0.031 


0.032 


0.021 


0.035 


0.037 


Total Error 


0.036 


0.035 


0.023 


0.039 


0.040 



Table 7: Systematic uncertainties in Am^ and in the mistag rates measured with semileptonic B 
decays for the likelihood fit. 



Source 


Arud 


Lepton 


Kaon 


NTl 


NT2 




[h ps-i] 










At Resolution 


0.012 


0.005 


0.009 


0.012 


0.005 


Background At 


0.002 


0.002 


0.002 


0.002 


0.002 


Background Resolution 


0.002 


0.002 


0.002 


0.002 


0.002 


Background Dilutions 


0.006 


0.008 


0.013 


0.026 


0.031 


Background Fractions 


0.006 


0.009 


0.011 


0.017 


0.032 


B'^ Backgrounds 


0.010 


0.009 


0.010 


0.004 


0.003 


B^ lifetime 


0.006 


0.001 


0.001 


0.001 


0.001 


z scale 


0.005 










z boost 


0.003 










Monte Carlo Correction 


+0.008 


-0.010 


-0.001 


-0.002 


-0.006 




± 0.009 


± 0.008 


± 0.006 


± 0.011 


± 0.011 


Total Systematic Error 


0.022 


0.018 


0.023 


0.035 


0.046 


Statistical Error 


0.020 


0.020 


0.016 


0.028 


0.025 


Total Error 


0.030 


0.027 


0.031 


0.045 


0.052 
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optimized At interval. Taking into account the BABAR vertex resolution, the optimum interval is 
found to be \At\ < 2.5 ps. This is so because the number of mixed events in this time interval 
is dominated by the mistag rate rather than by B^B^ mixing. Events with |At| > 2.5 ps have 
on average equal numbers of mixed and unmixed events due to B^B^ oscillations, and therefore 
contribute nothing to the determination of the mistag rate. We refer to this time- integrated method 
using a single optimized At interval as the "single-bin" method. 

The single-bin analysis uses the reconstructed hadronic and semileptonic B^ sample described 
in Section The number of tagged events in each category were summarized in Table ^ and 
Table ^. The background fractions were re-evaluated for the sample of tagged signal events with 
|At| < 2.5ps. 

To correct for the presence of backgrounds, we must add to Equation |6| a term to account for 
the contribution of each background source to the fraction of mixed events in the sample: 

Xohs = + (1 - 2xd) //3X/3, (16) 

P 

where fs^fp are the fraction of signal and each background source, respectively, xp is the fraction 
of mixed events in each background source, and Xobs is the observed fraction of mixed events. We 
restrict the sample to events with | At| < 2.5 ps; then Xd must be modified to represent the integrated 
mixing probability for | At| < 2.5 ps. Using the world-average value for Am^, 0.472zb0.017 U ps~^ Q, 
and taking into account the detector resolution function TZ{At), we find 

, 1. /-Z5;;{e-l^^l/^co.(Am,At) n{At)}d{At) ^ 
\d = - 1 — nrE nrrr, = 0.079. (17) 



Solving for Eq. 16 for w, and using the calculated value for x^, we obtain 



Xobs - fsX'd - J2i3 fpXp _s 



8.2 Results 

8.2.1 Hadronic sample 

The selection of hadronic B mesons was described in Section We use all tagged events with 
I At I < 2.5 ps, and determine the background fraction in the signal sample from a fit to the tit-es dis- 
tribution as described in Section |^. The signal region is defined as events with tties > 5.27GeV/c^. 
The fraction of mixed events in the background is determined by tag category using the sideband 
control sample, rriEs < 5.27GeV/c^. 

We use Eq. ^ to solve for the mistag rate for each tag category, obtaining the results shown in 
Table g 



8.2.2 Semileptonic sample 



The selection oi B ^ D*lv events was described in Section 6.3. We use tagged events with |At| < 
2.5 ps and re-evaluate the backgrounds for events in this time interval, with the control samples 
described above. The backgrounds are evaluated for each tag category and for each decay 
mode. The mistag fractions are calculated individually by tag category and decay mode using 
Eq. 18, and the results for the different decay modes are combined, using the statistical errors 
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Table 8: Mistag rate w and tagging separation Q as measured by the single-bin method in the 
hadronic and semileptonic B event samples. The mistag rates include small corrections corre- 
sponding to the difference between the generated and reconstructed values in simulated signal 
events. 



Tagging Category 


Hadronic 


Semileptonic 




Mistag Rate w Q = 


etag{l-2w)^ 


Mistag Rate w Q = 


etagil - 2wf 


Lepton 


0.120± 0.032 


0.061 


0.095± 0.018 


0.079 


Kaon 


0.207± 0.021 


0.126 


0.177± 0.014 


0.154 


NTl 


0.127± 0.035 


0.067 


0.200± 0.024 


0.041 


NT2 


0.342± 0.036 


0.016 


0.346± 0.024 


0.016 


All 




0.270 


— 0.290 



to weight the individual results. All systematic errors are conservatively taken to be correlated 
between the different decay modes. The determination of systematic errors will be discussed in the 
next section. The results are summarized in Table 

8.3 Systematic errors 

The various sources of systematic error in the single-bin method that have been investigated for the 
hadronic and semileptonic B samples are summarized in Tables ^ and |l^ and discussed in detail 
below. 

Table 9: Sources of systematic error for the mistag measurement on the hadronic sample in the 
single-bin method. See text for details. 



Type 


Variation 


Lepton 


Kaon 


MTl 


NT2 


Background 


±la 


0.004 


0.004 


0.008 


0.006 


Xd 


0.174 ±0.009 


0.005 


0.004 


0.005 


0.004 


Resolution 


see text 


0.002 


0.002 


0.002 


0.002 


wrong tag resolution 


see text 


0.005 


0.009 


0.007 


0.013 


Monte Carlo correction 




+0.009 


+0.004 


-0.026 


-0.007 


Monte Carlo statistics 




0.012 


0.008 


0.017 


0.015 


Total 




0.015 


0.014 


0.021 


0.021 



The systematic error due to background in the data samples is taken by varying both the 
background fractions and the fraction of mixed events associated to each background source, x/s- 
These quantities are varied by one standard deviation of the values measured in the background 
control samples described in Section ^. For the semileptonic sample, this is the dominant source of 
systematic error, primarily due to the limited statistics of the background control samples. 

The systematics uncertainties introduced by background from the decay B~^ — > D*~Xi~^i' are 
obtained by varying the fraction described in Section |6.3| as well as the mistag fraction of the 
charged B meson measured on data. 

The assumed At resolution function is a double-Gaussian where the initial parameters are taken 
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Table 10: Sources of systematic error for the mistag measurement from the semileptonic sample in 
the single-bin method. See text for details. 



Type 


Variation 


Lepton 


Kaon 


NTl 


NT2 


Background 


±la 


0.008 


0.010 


0.023 


0.020 


B~ D*-Xi+u 


see text 


0.009 


0.009 


0.009 


0.008 


Xd 


0.174 ±0.009 


0.012 


0.009 


0.009 


0.004 


Resolution 


see text 


0.002 


0.002 


0.002 


0.002 


wrong tag resolution 


see text 


0.005 


0.009 


0.009 


0.013 


Monte Carlo correction 




0.001 


-0.001 


-0.003 


-0.014 


Monte Carlo statistics 




0.010 


0.006 


0.013 


0.012 


Total 




0.021 


0.020 


0.031 


0.028 



from simulation. The two widths are multiplied by two scale factors which are determined in a fit 
to the data along with the fraction of events in the narrow Gaussian. The double- Gaussian RMS 
is consistent with that of the resolution function considered in Section |^ These three parameters 
are varied in a conservative way within the errors of the fit to the data to determine the systematic 
uncertainty on the mistag fraction, which is quite small. 

Finally, we consider the possibility that wrong tags have worse At resolution than correct tags. 
This effect has been studied in Monte Carlo simulation, where we observe a slightly larger RMS in 
the At distribution for events with wrong-sign tags. We assign a systematic error on the mistag rates 
by taking the default resolution function and changing the parameters of the resolution function 
for wrong tags by a corresponding amount. 

9 Comparison of likelihood and single-bin methods 

Combining the results obtained for the hadronic and semileptonic B samples for the likelihood fit 
method described in Section |^ and for the single-bin method described in Section ^ and taking into 
account the correlated systematic errors, we obtain the preliminary mistag rate results summarized 
in Table |ll[ The Monte Carlo corrections to the likelihood fit results summarized in Tables ^ and 
have been applied to the likelihood fit results of Table |5[ 

Table 11: Combined mistag results for the hadronic and semileptonic B samples, for the likelihood 
and single-bin methods. 



Tag Category 


Mists 

Likelihood 


Ig, w 

Single-Bin 


Lepton 
Kaon 
NTl 
NT2 


0.096 ± 0.017 ± 0.013 
0.197 ± 0.013 ± 0.011 
0.167 ± 0.022 ± 0.020 
0.331 ± 0.021 ± 0.021 


0.102± 0.016 ± 0.015 
0.187± 0.012 ± 0.015 
0.176± 0.020 ± 0.024 
0.345± 0.020 ± 0.023 



The single-bin fit results use a sub-sample of the events used in the likelihood fit, so the two 
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sets of results are correlated. The systematic errors differ, due to the different sensitivities of the 
two methods. Overall, the two methods agree well. The effective flavor tagging efficiency of the 
algorithm described in Section || is found to be Q w 0.28. 

10 Summary 

In 8.9 fb^^ of e^e^ annihilation data collected near the T{AS) resonance, we have obtained a prelim- 
inary measurement of the time-dependent B^B^ oscillation frequency using a sample of B^ mesons 
reconstructed in hadronic decay channels and in the semileptonic decay mode B^ D*^l'^u. 
From the hadronic B^ sample we measure the B^B^ oscillation frequency: 

Amd = 0.516 ± 0.031 (stat.) ± 0.018 (syst.) h ps~^ 

From the sample we measure the B^B^ oscillation frequency: 

Amrf = 0.508 ± 0.020 (stat.) ± 0.022 (syst.) h ps~^ 

Combining the Am^ results from the hadronic and semileptonic B samples, we obtain the prelim- 
inary result : 

Am^ = 0.512± 0.017 (stat.) ± 0.022 (syst.) fi ps"^ 

In combining the two results, we have taken all systematic error contribution to be fully correlated 
with the exception of the contribution due to Monte Carlo simulation statistics. 
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